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Abstract 

We report the discovery of a peculiar horizontal branch (HB) in NGC 6440 and NGC 6569, 
two massive and metal-rich Galactic globular clusters (GGCs) located in the Galactic bulge, 
within 4kpc from the Galactic Center. In both clusters, two distinct clumps are detected at 
the level of the cluster HB, separated by only ~ 0.1 magnitudes in the K s band. They were 
detected with IR photometric data collected with the "VISTA Variables in the Via Lactea" 
(VVV) Survey, and confirmed in independent IR catalogs available in the literature, and HST 
optical photometry. Our analysis demonstrates that these clumps are real cluster features, not 
a product of field contamination or interstellar reddening. The observed split HBs could be 
a signature of two stellar sub-populations with different chemical composition and/or age, as 
recently found in Terzan5, but it cannot be excluded that they are caused by evolutionary 
effects, in particular for NGC 6440. This interpretation, however, requires an anomalously high 
helium content (Y > 0.30). Our discovery suggests that such a peculiar HB morphology could 
be a common feature of massive, metal-rich bulge GGCs. 

Subject headings: globular clusters: general 



1. Introduction 

Our understanding of the complexity of Galac- 
tic Globular Clusters (GGCs) has impressively 
expanded in the last decade, propelled by the 
discovery that they can host multiple popula- 
tions of stars with a different c hemical enrich- 
ment history ( Piotto et al. 20051) . The classi- 
cal text-book definition of GGCs as prototypes 
of a simple stellar population, i.e. a chemically 
homogeneous aggregate of coeval stars, is now 
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out-dated. While a certain degree of inhomo- 
geneity of light c hemical elements is o bserved in 



nearly all GGCs (jCarretta et al. 2009), a spread 



a few very massive objects 
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; Cohen et al. 2010). 


Ferraro et al. 1 (2009) discovered two horizontal 



branches (HBs) in the Bulge GGC Terzan 5, sepa- 
rated by 0.3 magnitudes in the K s band. The exis- 
tence of multi-modality in the morphology of HBs 
has b een known for nearly four decades ( Harris! 



19751 ). and has been associated with the presence 
of mu ltiple stellar population s since shortly there- 
after (|Rood fc Crocker! 1 19851) . However, to date, 



Terzan 5 is the only GC known to have two dis- 
tinct HBs. The two features in Terzan 5 have 
a different spatial distribution, the brighter one 
being more centrally conce ntrated, more meta l 
rich ( A[Fe/H]~ 0.5 dex, lOriglia et all l201ll) 
and possibly helium enh anced (ID'Antona et al 

2009) 



2010T) and/or youn ger (jFerraro et al 



Lanzoni et al. I (l2010l) confirmed that Terzan 5 is 
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more massive than previously thought, and it 
could be the relic of a Bulge building block. 

In this Letter, we show evidence that the Bulge 
GGCs NGC 6440 and NGC 6569 host split HBs, 
similar to that of Terzan5. NGC 6440 is a high - 



metaUicity ([Fe/H]« -0.5. lOriglia et al. 1 120081 ) 
cluster, located 8.5 kpc from the Sun and onl 
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1.3 kpc from the Galactic center (jHarris 
2010 edition, H10). NGC 6569 is slightly less 



metal-rich ([Pe/H]w -0.79. IValenti et al. 1120111 ) 



and is found at a distance of 10.9 kpc from the Sun 
and 3.1 kpc from the Galactic center (H10). Both 
NGC 6440 and NGC 6569 are among the ten most 
luminous of the 64 GGCs located within 4 kpc 
from the Galactic center. 

2. Observations and reductions 



The "VISTA Variables in 
(VVV) Survey (|Minniti et al. 



t he V ia 
2010) 



Lactea" 
is one of 

the six ESO Public Surveys operating on the 4- 
meter Visible and Infrared Survey Telescope for 
Astronomy (VISTA). VVV is scanning the Galac- 
tic bulge and the adjacent part of the southern 
disk (-65 < / < -10, -2 < b < +2), in five near- 
IR bands (YZJHK S ) with the VIRCAM cam- 
era 



( Emerson et al. 201dl ). an array of sixteen 



2048x2048 pixel detectors with a pixel scale of 
O'.'341/pix. VVV images extend four magnitudes 
deeper and exhibi t increased spatial resolution 
(jSaito et al. l l2010h versus Two M icron All Sky 
Survey (2MASS. ISkrutskie et al~1 120061) . which 



is particularly important for mitigating contam- 
inated photometry in crowded regions near the 
Galactic center. 

We retrieved from the Vista Science Archive 
website@ the VVV images of the two GGCs, 
pre-reduced at the Cambridge Astronomical Sur- 
yey Unit (CASUfl with the VIRCAM pipeline 
(Irwin et al. 2004). The selected data consist 



of four frames, sampling twice each point in an 
area of 17' x 22' around the GGCs, in each of 
the ZYJHK S filters, plus 17 and 11 additi onal 
epochs in the K s passband ( Saito et al. 1120121) for 
NGC 6440 and NGC 6569, respectively. The VVV 
images of the two clusters, extracted from a single 
K s frame, are shown in Figure [TJ 
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The PSF-fitting photometry was obtained with 
the VVV-SkZ_pipeline (VSp, Mauro et al. sub- 
mitted), code based on DAO PHOT and ALL- 
FRAME (jStetson I [19871 Il994l ) procedures, opti- 
mized for the VVV data. The photometry was 
tie d to 2MASS JHK s stan dards, as described 
in Moni Bidin et al. I (|201lh and IChene et al. 
(|2012h . Combining all the 36 and 24 K s mea- 
surements, the final photometric errors were 0.003 
and 0.008 mag respectively at the brightness level 
of the cluster HB. The K s errors for NGC 6569 
are costant with distance, while in NGC 6440 they 
increase up to 0.005 mag for r < 0.'7. The com- 
pleteness of our photometry is heavily affected 
by crowding in the inner 0'.7 of both clusters, as 
can be appreciated in Figure Q] and [4] Due to 
incompleteness, ~ 80% of the detected HB stars 
stay outside this problematic inner region, where 
crowding is not a significant issue. 

3. Results 

The (J - K S ;K S ) Hess diagrams (HD) for the 
stars detected within 1.'83 and 1.'66 from the clus- 
ter center, respectively, are shown in Figure [2] and 
[3l They were obtained calculating the number 
of stars in a bin of width 0.06 mag in J — K s 
(0.04 mag in the lower panels) and 0.04 mag in 
K Sl moved along the axis with steps of 0.015 mag 
in color (0.01 mag in the lower diagram) and 
0.01 mag in magnitude. The HDs of both GGCs 
reveal a peculiar HB morphology, with two dis- 
tinct clumps separated by ~ 0.1 mag. They will 
be referred to as HB-A and HB-B (lower panels of 
Figures [2] and [3|) for the brighter and fainter one, 
respectively. The overdensity observed at redder 
color is the red giant branch (RGB) bump: in fact, 
Valenti et al~l (|2005l V05) found K Stbump =U.08 
in NGC 6569, and their K S: b um p-{Fe/H] relation 
predicts K s ^ ump «14.1 in NGC 6440, in good 



agreement with our data. 

To verify if both HBs are real and belong to 
their host cluster, we checked the data for stochas- 
tic fluctuations as a cause of the overdensity, and 
analyzed their spatial distribution as a function of 
the central distance. Furthermore, we compared 



our da ta with the IR p hotometry of lValenti et al. 



(|2004L V0 41 andlVOa. and w ith the optical HST 
data from lPiotto et all (|2002h . 
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Fig. 1. — 2' x 2' linear gray-scale if s -band VVV 
images of NGC6440 (upper panel) and NGC6569 
(lower panel). The circle indicates the 0'7 inner 
limit of photometric completeness commented in 
the text. 



Dereddening. We used the maps from lGonzalez et al. 
( 201l[ ) to correct for reddening. They reveal that 
E(J - K s ) - 0.5 - 0.7 in the r = 1.'8 field of 
NGC6440 under analysis. The case of NGC6569 
is much less extreme, with E(J—K S ) = 0.20—0.24. 
The HD of NGC 6440 shows a clear improvement 
(see Figure [5]) , with the two features less blurred 
and HB-B still presenting a slope. For NGC 6569 
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Fig. 2. — (J — K s ; K s ) Hess diagrams of HB area 
of NGC 6440. The dereddened version is presented 
in the insert plot in upper panel. In lower panel, 
the groups HB-A and HB-B discussed in the text 
are spotlighted. The arrows show the effect of the 
reddening. 



the dereddened HD is approximately similar to 
the raw one, as expected. 

Checking for stochastic variation. We reran 
the procedure on four subsets of the original data, 
each one containing the ZY and JHK S data, but 
different K s epochs: one subset included only the 
first epoch, while a unique set of three epochs were 
used in each of the three following subsets. The 
declared photometric errors in K s passband vary 
from 0.007-0.009 mag to 0.003-0.005 mag at the 
level of the HB. Comparing the (J — K s ; K s ) HDs, 
obtained with the previous spacial selection and 
sampling procedure, both GGCs always exhibit a 
split HB, with only negligible differences in their 
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of 2" for NGC 6440, while for NGC 6569 we used 
the values of 15" and 6", respectively. The SDs of 
the two groups steeply decay at increasing radii, in 
both GGCs, and their members are distributed on 
the CCD with circular symmetry. The radial pro- 
file of the two features in NGC 6569 is identical. 
The HB-B group in NGC 6440 is more populated 
than the brighter HB-A by a factor of two, but a 
Kolmogorov-Smirnov test reveals that their radial 
behavior coincides also in this case. The stellar 
counts in the inner 0.'7 ~ 5r c of NGC 6440 (where 
r c is the core radius from H10) are incomplete be- 
cause of crowding. The photometry of NGC 6569 
is also incomplete for r < 0'7 ~ 2r c . 

Th e radia l profile of NGC 6569 was fit with a 



King I (|1962f) profile of the form 



i+ - 



i+ - 



-F. 



(1) 

where k is a scale parameter, r t is the tidal radius 
and F the field contribution. For NGC 6440, we 
used the approximation for r r c 



Fig. 3.— Same as Figure[2]for NGC 6569. 



morphology. As an additional test, we checked 
the HBs of other GGCs in the VVV, namely 
NGC 6380, NGC 6441, NGC 6528 and NGC 6553, 
finding no evidence of a split or peculiar HB. 



f(r) = krl 



F . 



(2) 



For NGC 6440, the fit leads to r t ,A = 5fl ± Q'.7 and 
Tt,B — 5.'2 ± 0.'5, consistent with r t — 5'84 quoted 
by H10. The core radius and the scale parameter 
cannot be separated and estimated individually. 
The SD of the two features in NGC 6569 are com- 
patible with r c = 0f35 and r t = 7.15 (H10). 



Field Contamination and Spatial Distribu- 
tion. We checked the field contribution to the 
HDs, selecting an annular region with the same 
area of the previous selection, but just outside the 
tidal radius. The HDs of the field are barely pop- 
ulated at the HB location, and the field contami- 
nation is negligible. This results is evident even in 
Figure HI where the number counts drop to near- 
zero levels at large distances from the cluster cen- 
ters. 

The behavior of the stellar densities (SDs) with 
distance r from the center is shown in Figure 2] for 
the two features highlighted in the lower panels of 
Figures[2]and[3] The SDs (stars per arcmin 2 ) were 
calculated with a bin width of 10" moved at steps 



Comparison with previous photometry. 

We matched our VVV photometry of NGC 6440 
and NGC 6569 with the catalogs of IV04j a nd 
V05l . respectively. The photometry of |V04| is 
based on observations with the near-IR camera 
IRAC2@ESO/MPI 2.2m, covering a 250" x 250" 
field center ed on the cluster. Similarly, the pho- 
tometry of was performed on data collected 
with the near-IR camera SOFI@ESO/NTT, , cov- 
ering a 300" x 300" field centered on the clus- 
ter. The estimated internal photometric errors 
are lower than 0.03 mag. Both photometries were 
calibrated onto the 2MASS photometric system 
and astrometrically corrected by using the 2MASS 
catalog. 
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-?^s,V05 = 14.26 and i4T S! v 5 = 14.36, respectively, 
for the HB-A and HB-B clumps in NGC 6569, with 
a dispersion of a = 0.07, and /sT s yy\/=14.26 and 
14.35 for the same features in our photometry, 
with a dispersion of a = 0.02. Thus, the mean 
magnitude of the clumps of both GGCs is identi- 
cal in VVV and Valenti et al.'s catalogs, but the 
separation is four times more statistically signif- 
icant in the VVV data. This result proves that 
the HBs of both GGCs are intrinsically split in 
magnitude, with a bright er and a f ainter part that 
remain separated even in IV04l and IV05l photome- 
try, respectively, once the stars are identified. 

The strong differential reddening affecting the 
region of NGC 6440 causes the HB to be strongly 
tilted at optical wavelengths, hence a simple lumi- 
nosity distribution does not show a bimodal be- 
havior. For t his reason, we ana lyzed the HST 
optical data of IPiotto et al. ( 20021 ) projecting the 
position of each HB star along the HB slope, ac- 
cording to the equation 



0.2 0.4 0.6 



2.2 2.4 2.6 



F555J< r =F555W„ r - a[(.F439W^ 

-(F439M/„ r -F555W„ r ) ], 



F555W nr ) 
(3) 



Fig. 4. — The distributions of the stellar density 
as function of the distance (bin width of 10"moved 
with steps of 2") for the HB-A (solid line) and HB- 
B (dashed line) features. The error bars are plot- 
ted only every 0.'3 and every 0.'25 for NGC 6440 
and NGC 6569, respectively; the King profile fits 
are overplotted. The vertical lines indicate the 
limit for the incompleteness. 



For both GGCs, the lumi nosity dis tribut ions 
in the K s mag nitudes of the Ivoi and [vol cat 



alogs for the matched stars do not show a clear 
bimodal distribution. However, when the stars be- 
longing to the HB-A and HB-B groups are iden- 
tified, their luminosity distributions are different, 
as shown in Figure [5] For NGC 6440, the Gaus- 
sian fits of the two distributions are centered at 
K s yo4 = 13.55 and if s ,vo4 = 13.66 for HB- 
A and HB-B, respectively, with a dispersion of 
a = 0.12, while in our VVV photometry the values 
are -ftT s yyv/=13.55 and 13.67, respectively, with 
a dispersion of a = 0.03. Analogously, we find 



where a — 3.7 is the slope of the HB and 
(F439W„ r - F555W nr ) = 2.2 is the HB mean 
color. As advised by the authors, the magnitudes 
adopted were those not corrected for reddening 
(nr). To avoid contamination from RGB stars, we 
selected only the sources with 18.2 < FhhhW nr < 
19.1 and 2 < (_F439W„ r - F555W nr ) < 2.4. 
The distribution of ^5551^,,, calculated with bin 
width of 0.08 mag and step of 0.04 mag (see Fig- 
ure E|) , reveals a clear double peak separated by 
~ 0.23 mag, with the fainter peak 1.6-1.7 times 
more populated. A similar procedure was per- 
formed for NGC 6569 also, but we were not able 
to disentangle any bimodality. 

4. Discussions and Conclusions 

The analysis of VVV data reveals that the HB 
of the GGCs NGC 6440 and NGC 6569 is split into 
two distinct clumps. This behavior is not intro- 
duced by stochastic fluctuations of the density in 
the CMD, or induced by photometric errors, as it 
is found to be identical in four independent sub- 
sets of data. Field contamination is not the cause 
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Fig. 5. — Luminosity distribution in K s vo4 and 
K s vo5 passband (in upper and middle frame for 
NGC6440 and NGC 6569, respectively) for the 
matched stars belonging to the two features HB- 
A and HB-B of the respective GGCs. Luminosity 
distribution of the values F555W£ r for NGC 6440 
in the lower frame. The center and the sigma re- 
sulting from the Gaussian fits are also shown. 



either, because the members of both the HBs are 
distributed with spherical symmetry with respect 
to the cluster center, and their density steeply de- 
cays with distance. The separation in NGC 6440 is 
even cleaner after applying a differential reddening 
correction, while in NGC 6569 it remains similar, 
presenting lower differential reddening. This HB 
split is, however, not found in the same VVV data 
of four less massive Bulge GGCs. 

The magnitude difference between the two HB 
clumps is only ~ 0.08 — 0.1 magnitudes in K Sl 
smaller than in Terzan5 by a factor of three. It 
is thus not surprising that this HB split passed 
unnoticed in previous investigations, also con- 
sidering the strong differential reddening affect- 
ing the NGC 6440 fie ld (AE(J - K s ) = 0.2, 



iGonzalez et al. 2011 ). We find that the di 



chotomy is blurred by obse rvati onal errors in the 
IR photometry of|V0J and IV05l . but the two fea- 
tures are clearly separated even in their data, once 
their members are identified in their catalog. The 
HST optical data of NGC 6440 from IPiotto et al 



(|2002l ) show two peaks separated by 0.23 mag 
in the corrected magnitude F555W„ r defined in 
Equation [3J 

The fainter HB of NGC 6440 is bluer than 
the brighter one. This rese mbles what was pre- 



viously found in Terzan5 (|Ferraro et al. I [2009; 



Lanzoni et al. 2010l ). In addition, the fainter HB 



of NGC 6440 is more populated than the brighter 
one by about a factor of two, slightly higher 
than what was found in the centra l regions of 



Terzan5 (~ 1.6, Ferraro et al. 20091 ). However 



these results are not directly comparable, because 
our photometry is incomplete in the inner 0.'7 of 
NGC 6440. HST data suggest that in the central 
region this ratio could be lower (~ 1.6), as ex- 
pected if, analogously to Terzan 5, the brighter 
HB was more centrally concentrated. On the other 
hand, we did not detect any difference in the ra- 
dial density profile of the two clumps, so the issue 
remains open. 

The two HBs of Terzan 5 are associated with 
two sub-populations of different mctallicity, with 
the brighter HB being richer in iron by ~ 0.5 
dex (iFerraro et al. ll2009HOriglia et al.| [201lh. Ac- 



cording to ISalaris fc Girardi (l2002h . a difference 
of ^0.3 dex would be expected, if the observed 
HB splits are interpreted only in terms of metal- 
licity. Nevertheless, this is only a rough upper 
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limit, because differences in helium content and 
age also can contr ibute to cause the same split. 
Origlia et al. ( 20081 ) measured the metallicity of 
ten stars in NGC 6440, finding a dispersion of only 
0.06 dex, compatible with observational errors. 
However, their targets are mainly located within 
0.'7 of the center, where the HB-A members could 
be few if the population ratio is constant at all ra- 
dial distances. Hence, Origlia et al.'s sample likely 
contains only a small quantity of HB-A stars, and 
their results are insufficient to exclude a metallic- 
ity spread in this cluster. 

Contrary to the case of NGC 6440 and Terzan5, 
the two groups identified in the HB of NGC 6569 
have approximately the same color and the same 
radial distribution, the fainter HB B being 1.3 
times more populat ed than the other clump. 
Valenti et al. I ( 2011 ) measured the metallicity of 



components are well separated and with a narrow 
color s pread, is very different to the simu lated 



six stars in this cluster, finding a bimodal distri- 
bution with two groups separated by ~ 0.08 dex. 

It is possible that the peculiar HB morphol- 
ogy discovered in NGC 6440 is a pure evolution- 
ary effect, and not a signature of the presence 
of sub-populations. In fact, wedge-shaped HBs 
are predicted und er special circumstances, as de- 



picted in Fig. 4 of ICatelan fe de Freitas Pacheco . 
(|l996l) and Fig. 1 of iDorman et al.l (|l989t ). and 



statistical effects could lead to the actual bimodal 
distribution of HB magnitudes. These features 
are found in the luminosity-temperature plane, 
but the simulated optical CMDs reveal only a 
clump at the bluer (and brighter) end of the 
sloped HB, at variance with what is observed 
in HST data, while the behavior of these fea- 
tures in the IR bands has not been simulated. 
Hence, this interpretation seems unlikely, but it 
cannot be excluded and represents an intriguing 
possibility. In fact, the high metallicity alone can- 
not explain the formation of a wedge-like HB, 
and a very high helium abundance (Y > 0.30) 
is required. Such a He-enriched field population 
has been recently su ggested in the bulge (e.g., 
Nataf fe Gould 1 120121 ). although at higher metal- 
licities. Hence, interpreting the HB morphology 
of NGC 6440 as an evolutionary effect implies that 
the helium content of this cluster must be anoma- 
lously high, actually higher than what is predicted 
at [Fc/H]=— 0.5 b y the models of Bulge chemical 
enric hment (e.g., ICatelan fc de Freitas Pacheco 
1996f ). The split HB observed in NGC 6569, whose 



color s pread, is very dittcrcnt to tne simulated 
HBs of ICatelan fc de Freitas Pacheco" (1996) and 
IDorman et all (|l989h . The interpretation of these 
features as an evolutionary effect induced by a 
high helium content is unlikely 

The fainter HB of NGC 6440 (HB-B) is tilted, 
the brightness of its stars increasing at bluer col- 
ors. This is clearly visible in Figure [5J where 
we indicate the direction of the redden ing in the 
IR bands from ICatelan et. al (l201ll) for com- 
parison. Dereddening the photometry, the slope 
is still present, but the map resolution of 1' 
does not permit strong claims. The tilt is more 
pronounced i n the optical HST photometry of 
IPiotto et all (|2002t ). where we measured a slope 
A(V)/A(B - V) ss3.7, which is higher than the 
standard reddening law Ry = Ay / E(B — V) «3.1. 
Moreover, the optical extinction is non-standard 
toward the G alactic bulge, and Ry can be as 



low as ~2.5 ( Nataf et al. I 120121) . In conclu- 



sion, the slope of the HB in NGC 6440 is di- 
rected approximately aligned with the redden- 
ing vector, but it is steeper than the expecta- 
tions of interstellar reddening. This behavior 
was already observed in NGC 6388 and NGC 6441 
(jSweigart fc Catelarj[l995 iBusso et al.ll2007t ). two 
other massive, metal-rich Bulge GGCs, and it 
was attr ibuted to an anomalou sly high helium 
content ( Caloi fc D 'Antonall2007h . Very interest- 
ingly, the HB stars of NGC 6388 could show the 
same peculiar properties observed in u Centauri 
(|Moehler fc SweigartteOoHlMoni Bidin et all 20 HI ), 
the most famous cluster hosting a complex mix of 
sub-populations with different chemical enrich- 
ment histories. 

Our results indicate that Terzan 5 is not a 
unique object. A complex HB morphology could 
be a relatively common feature among metal-rich, 
massive Bulge GGCs. This is not detected in less 
massive objects (e.g. NGC 6528, NGC 6553), nor 
in equally massive but metal-poor Bulge GGCs, 
such as M 22 and M 28, whose HB is very extended 
toward the blue. T he large metallicity spread 



observed in Ter zan 5 ( Origlia et al 



present in M22 (|Marino et all 12009) . Further in- 



20111) is also 



vestigations are needed to unveil if the HB splits 
reflect the presence of two stellar populations with 
different chemical composition and/or age. 
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